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Abstract 

An alternate Hamiltonian H different from Ostrogradski's one is 
found for the Lagrangian L = L{q,q,q), where d'^L/d{q)^ 7^ 0. We 
add a suitable divergence to L and insert a = q sen. d b = q. Contrary 
to other approaches no constraint is needed because a = his one of the 
canonical equations. Another canonical equation becomes equivalent 
to the fourth-order Euler-Lagrange equation of L. Usually, H be- 
comes quadratic in the momenta, whereas the Ostrogradski approach 
has Hamiltonians always linear in the momenta. 

For non-linear L = F{R), G = dF/dR 7^ the Lagrangians L and 
L = F{R) with F = 2R/G^ - ?>L/G^, gij = gij and R = 3R/G^ - 
AL/G^ give conformally equivalent fourth-order field equations being 
dual to each other. This generalizes Buchdahl's result for L = R?. 

The exact fourth-order gravity cosmological solutions found by 
Accioly and Chimento are interpreted from the viewpoint of the in- 
stability of fourth-order theories and how they transform under this 
duality. 

Finally, the alternate Hamiltonian is applied to deduce the Whee- 
ler-De Witt equation for fourth-order gravity models more systema- 
tically than before. 
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1 Introduction 

Higher-order theories, especially fourth-order gravity theories, are subject 
to conflicting facts: On the one hand, they appear quite naturally from 
generally accepted principles; on the other hand, they are unstable and so, 
they should be considered unphysical. 

Here, we want to attack this conflict from two directions. First: The 
Ostrogradski approach [1] to find a Hamiltonian formulation for a higher- 
order theory is the most famous (see e.g. refs. [1-9]) but possibly not the 
best method. To check this hypothesis, we present an alternate Hamiltonian 
formalism for fourth-order theories in set. 2. It systematizes what has been 
sporadically done in the literature for special examples. 

Set. 3 deals with fourth-order gravity following from a non-linear La- 
grangian L{R). The conformal equivalence of these theories to theories of 
other types is widely known, but the conformal equivalence of these theories 
to theories of the same type but essentially different Lagrangian is much less 
known. We fill this gap by proving a duality theorem between pairs of such 
fourth-order theories in subsection 3.1. The instabihty of these theories from 
the point of view of the Cauchy problem is subject of subsection 3.2. 

Set. 2 applies to arbitrary theories, set. 3 to gravity, and both are applied 
to fourth-order cosmology in set. 4. We re-interpret known exact solutions 
(Priedmann models in subsection 4.1 and Kantowski-Sachs models in 4.2) 
under the stability criteria mentioned before. 

In the final set. 5 we discuss quantum effects and give hints (which shall 
be outlined in a future paper with S. Renter) how to apply the alternate 
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Hamiltonian formalism of set. 2 to the Wheeler-De Witt equation for a 
cosmological minisuperspace model within fourth-order gravity. 



- ^ („)^ („) (1-1) 



The rest of this introduction shortly reviews papers on higher-order 
theories. Eliezer and Woodard [1] and Jaen, Llosa and Molina [2] represent 
standard papers for the generalization of the Ostrogradski approach to non- 
local systems (see also [3]) and to systems with constraints (see also [4 - 6]) 
applying Dirac's approach. 

Let the Lagrangian L be a function of the vector and its first n temporal 
derivatives Qa, Qa, ■ ■ ■ , q^^- The Hessian is 

dqa^dqp' 

and the non-vanishing of its determinant defines the regularity of L. In the 
following we do not write the subscript a; one can think of q as being a 
point particle in a (one- or higher-dimensional) space. In the Ostrogradski 
approach, Q = q is taken as additional position variable. This leads to an 
ambivalence of the procedure, because it is not trivial to see at which places 
q has to be replaced with Q, cf. [7]. We prevent this ambivalence in our 
alternate Hamiltonian, cf. set. 2, by putting Q = q. 

Ref. [8] discusses higher-order field theories. The problem is the lack of an 
energy bound, typically two kinds of oscillators with different signs of energy 
exist. Usually, one restricts the space of initial conditions to prevent negative 
energy solutions. The authors of ref. [8] redefine the energy analogous to 
the Timoshenko model, so one gets a positive mechanical energy inspite of 
an indefinite Ostrogradski Hamiltonian, they write: "An appealing aspect of 
this approach is the absence of any constraint." So it has this property in 
common with our approach set. 2, but it is otherwise a different one. 

A second standard procedure [2, 8, 9] for dealing with higher-order 
Lagrangians is to consider them as a sequence in a parameter e, so one can 
break the Euler-Lagrange-equation into a sequence of second order ones. In 
[9] this is called "reduction of higher-order Lagrangians by a formal power 
series in an ordering parameter." [9] deals also with the Lie-Konigs theo- 
rem: a local Hamiltonian is always possible, and they consider some global 
questions. 
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Let us repeat the famous counter-example [10] Douglas (1941): it is an 
example of a second order system not following from a Lagrangian: 

X + ij = Q y + y + ex = Q 

It follows from the Lagrangian 

L = - y2 + e{xy -xy- x^)\ 

for e 7^ and has no Lagrangian otherwise. We mention this example to 
show that the following recipe need to to work always. Recipe for higher- 
order theories: "Write down the Euler-Lagrange equations, break them into 
a sequence of second order ones by introducing further coordinates. Find 
Lagrangians for these second order equations." 

A powerful method for dealing with a classical Lagrangian 

L^\gijctq' -V{q) (1.2) 

is given in [11] and shall be applied in cosmological minisuperspace models 
hke in set. 5. The Euler-Lagrange equation to Lagrangian (1.2) reads 

and is fulfilled for geodesies in the Jacobi-metric 

9ij = {E-V)gij 

Remark: For constant potentials V this is trivial, for non-constant potentials 
the constant E must be correctly chosen to get the result, for E = V ii breaks 
down, of course. 

Stelle [12] cites Ostrogradski [1] but uses other methods to extract 
different spin modes for fourth-order gravity. In [13], a regular reduction 
of fourth-order gravity similar to the method with an ordering parameter 
mentioned above has been proposed as follows: In the Newtonian limit one 
has 

A$ + /3AA$ = 47rGp, 
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then one restricts to solutions which can be expanded into powers of the 
coupling parameter (5. Argument: If /3 is a parameter, this is well justified, 
if it is a universal constant, then this restriction is less satisfying. Comment: 
This restriction excludes the usual Yukawa-like potential ^exp(— r/y^), so 
one may doubt whether this method gives the right solutions. Let us further 
mention ref. [14] for non-local gravitational Lagrangians like L — i?n~^i?+A 
in two dimensions and refs. [15, 16] for the linearized i?^-theory. 

To facilitate the reading of set. 2, we pick up the example cq. (5) of [5]: 

L^lif + Aqe + Aq^y" (1.3) 

The equation of motion is [5, eq. (6)] 

2?^^) + I2qq^^'^ + 9{qf + 18^^^' = (1.4) 

A good check of the validity of the formalism is the following: For a constant 
c > and ^ > 0, each solution of 

q = -2q^ + c^q 

is also a solution of eq. (1.4). 

By adding a divergence to eq. (1.3) one gets L = (g)^e^'^. The alternate 
formalism requires to use q^ — q and q^ — q as new coordinates. So we get 

L=(g2)2 exp(3g^) (1.5) 

Eq. (1.5) represents the ultralocal Lagrangian mentioned in [5]. It is cor- 
rectly stated in [5] , that the alternate formalism does not work for this version 
eq. (1.5) of the system. This clarifies that the addition of a divergence to a 
higher-order Lagrangian sometimes influences the applicability of the alter- 
nate Hamiltonian formalism. So one should add a "suitable" total derivative 
to the Lagangian. "Suitable" means, that the space of solutions is the same 
at both sides, and that the relation between the various coordinates is en- 
sured without imposing any constraints. It turns out, that the Lagrangian 
L differing from L, eq. (1.3) by a divergence only 

L = -[(f + Q>qq^ + 2qq^^^]e^'^ (1.6) 
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does the job. Of course, the variations of L, L, and L with respect to q all give 
the same equation of motion (1.4). But only in version (1.6) the alternate 
formalism (insertion of the equation Q = q and then apply the usual formulas 
of classical mechanics - to avoid ambiguities with the square-sign we have 
replaced q^ by q and q^ by Q) leads correctly to the Hamiltonian [5, eq. (7)]: 

H = -^{pP -3P^Q)e-'^'' + Q^e^'' (1.7) 

It essentially differs from the Ostrogradski approach because terms only linear 
in the momenta do not appear; so one of the criteria for unboundedness of 
energy fails to be fulfilled. The integrability condition Q = q and the equation 
of motion (1.4) both follow from the canonical equations of eq. (1.7); no 
constraint is necessary to get this. 



2 The alternate Hamiltonian formalism 

Let us consider the Lagrangian 

L = Liq,q,q) (2.1) 
for a point particle q{t), a dot denoting ^ and 

The corresponding Euler-Lagrange equation reads 

dq dt dq dt"^ dq 
We suppose this Lagrangian to be non-degenerated, i.e., L is non-linear in q. 
The highest-order term of eq. (2.2) is 

therefore, non-degeneracy (= regularity, cf. eq. (1.1)) is equivalent to require 
that eq. (2.2) is of fourth order, i.e. 
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(If q is a vector consisting of m real components then this condition is to be 
written as Hessian determinant.) 

If we add the divergence ^G{q, q) to L, we do not alter the Euler- 
Lagrange equation (2.2). Furthermore, the expression -^G is linear in g, 
and so its addition to L does not influence the condition of non-dcgcneracy. 
The addition of such a divergence can therefore simply absorbed by a suitable 
redefinition of L. 

In the next two subsections we add a special and a more general diver- 
gence to get a Hamiltonian formulation different from Ostrogradski's one. In 
the preprint by Kaspcr [4] a similar consideration has been made at the La- 
grangian's level. Subsection 2.1 represents only a special case of subsection 
2.2, but we write it down, because it has the advantage that the formulas 
can be given explicitly, and so the formalism becomes more transparent. 



2.1 A special divergence 

The addition of the following divergence is no more done by a redefinition of 
L 

Ldiv^f^imqq], m^^O (2.3) 

and we consider L = L + L^iv The Euler-Lagrange equation is again eq. 
(2.2). Using 

we get 

L^L + f{q)eq + /(?)[(g)^ + qq^'^ (2.4) 

which contains third derivatives of q. 
We introduce new coordinates 

a — q , b — q (2.5) 

(In the Ostrogradski approach, the second coordinate is q, instead.) It is 
obvious that there is exactly this one compatibility condition: 

a = 6 (2.6) 
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Let us insert eq. (2.5) into eq. (2.4). This insertion becomes unique by 
the additional requirement that L does not depend on second and higher 
derivatives of a and 6, i.e., 

L — L{a, a, b, b) 

giving 

L = L{a, a, b) + f\a)a% + /(a) [b^ + ab] (2.7) 

(In the Ostrogradski approach, there remains an ambivalence which of the 
q in the original Lagrangian is to be interpreted as second coordinate and 
which as time derivative of the first one.) 

The momenta are defined as in classical mechanics by 

dL dL 

= fid ■ = ^ 

(In the Ostrogradski approach, an additional term is necessary.) Inserting 
eq. (2.7) into eqs. (2.8) we get 

dL 

Pa = -^ + 2f'{a)ab + f{a)b (2.9) 

and 

Pb = f{a)d (2.10) 
Because of /(a) ^ 0, cf. eq. (2.3), we can invert eq. (2.10) to 

Pb 



a — 



Inserting eq. (2.11) into eq. (2.9) and dividing by /(a) we get 



(2.11) 



^ >a-§-2/'(a)6^] (2.12) 



It is instructive to make a more general consideration: The question, whether 
eqs. (2.9, 10) can be inverted to d, 6, can be answered by calculating the 
Jacobian 

J ^ d{pa,pb) _ dpa dpb dpa dpb ^2 13) 

d{a,b) dd db db dd ^ ' ^ 

We insert eqs. (2.9, 10) into eq. (2.13) and get 

J = - [f{a)r (2.14) 
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Because of / 7^ one has also J 7^ and the inversion is possible. This more 
general consideration gave the additional information that the Jacobian is 
always negative; this may be the hint to some kind of instability. 
We define the Hamiltonian H as usual by 

H = apa + bpb - L 

i.e., with eq. (2.7) we get 

H^apa + bpb-L- f{a)a% - f{a) [6^ + ab] (2.15) 

Here we insert a according to eq. (2.11) and get the Hamiltonian H — 
H{a,pa,b,pb)- The factor of b in H automatically vanishes, so we do not 
need eq. (2.12). The canonical equations read 



further 



and 



and 



dH 
dpa 

dH 

dpb 



dH 
da 



= a (2.16) 

= b (2.17) 

= -Pa (2.18) 

= -Pb (2.19) 



dH 
'db 

The whole procedure is intended to give the following results: The Hamil- 
tonian H shall be considered to be a usual Hamiltonian for two interact- 
ing point particles a{t) and b{t). One of the canonical equations shall be 
equivalent to the compatibility condition eq. (2.6) and another one shall be 
equivalent to the original Euler-Lagrange equation (2.2), whereas the two 
remaining canonical equations are used to eliminate the momenta Pa and pb 
from the system. The next step is to find those Lagrangians L which make 
this procedure work. From eqs. (2.15) and (2.11) we get 

rr PaPb J I Pb Plf'{o)b .,2 /oon^ 
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In this form, eq. (2.16) coincides with eq. (2.11) and (2.17) with (2.12). So 
we may use eqs. (2.9, 10) in the following, because they are equivalent to 
eqs. (2.11, 12). 

Now, we use eqs. (2.19), cancel p& by use of eq. (2.10) and get 

= — + 26/(a)-a/(a) (2.21) 

In order that the compatibility relation eq. (2.6) follows automatically from 
eq. (2.21), one has to ensure that /(a) ^ (which is already assumed) and 
that 

o = ^ + 6/w 

identically takes place. The condition of non-degeneracy. 



is then also automatically fulfilled. One has the following possible Lagrangian 

L^~j{a)\?^K{a,a) (2.22) 

where K is an arbitrary function, but, for simplicity, we put K — 

The last of the four canonical equations to be used is eq. (2.18) reading 
now with eqs. (2.9, 10, 20) 

= /6 + 2/' hh + \f }? + /" a% (2.23) 

If we insert here eq. (2.5) we get exactly the same as the Euler-Lagrange 
equation (2.2) following from the Lagrangian 

(2.24) 

Result: For every Lagrangian of type (2.1) which can be brought into 
type (2.24) with / 7^ the addition of the divergence (2.3) makes it possible 
to apply the new coordinates (2.5). Then the system becomes equivalent to 
a classical Hamiltonian of two particles, and the relation (2.6) between them 
follows without imposing an additional constraint. 
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2.2 A general divergence 

In this subsection we try to generalize the result of the previous subsection 
by avoiding to prescribe the special structure (2.3) of the divergence to be 
added. We substitute eq. (2.3) by 

Ldiv = ^HQ: <i: q) (2-25) 

Keeping eqs. (2.5) we get instead of eq. (2.7) now 

L{a,a,b) + hid + h2b + h3b (2.26) 

where hn denotes the partial derivative of h with respect to its nth argument. 
Using eqs. (2.8), (2.10) is now replaced with 

P6 = /i3(a,d,6) (2.27) 

Eq. (2.13) is kept, and (2.14) is replaced with 

J=-(M' (2-28) 

We have to require that /i23 7^ 0, and then the equation — hz is locally 
invertible as a — F(pb,a,b). Prom this definition one immediately gets the 
identity Fi /i23 = 1- Two further identities to be used later are not so trivial 
to guess. To derive them, let us for a moment fix and then calculate the 
increase of hs with increasing a and b resp. The assumed constancy of 
yields the equations 

hi3 + F2h23 = (2.29) 

and 

h3 + F3/123 = (2.30) 

resp. to be used for deducing the generalization of eq. (2.21). One gets the 
result: For /i23 7^ (which is already presumed), the compatibility relation 
(2.6) follows automatically from the canonical equation (2.19) if and only if 

= La + /i2 (2.31) 

is identically fulfilled. One can see: The condition of non-degeneracy of the 
Lagrangian (2.1) namely 

^337^0 
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is equivalent to the condition /i23 7^ 0. For any given non-degenerate La- 
grangian we can find the appropriate divergence by solving eq. (2.31) as 
follows 

KlAA)^- L3{q,x,q)dx (2.32) 
Jo 

All other things are fully analogous: 

H^\pa- hi{a, F, b)]F - h2{a, F, h)h - L{a, F, h) (2.33) 

where F = F{pb,a,b). Eq. (2.19) with (2.30) gives the compatibility con- 
dition (2.6). Eq. (2.18) with (2.29) is equivalent to the Euler-Lagrange 
equation (2.2). 

Let us summarize this section: For the Lagrangian L = L{q, q, q) where 
d'^L/d{qY 7^ we define L = L + L^i^ where 




We insert a — q and h — q, define the momenta Pa = ^ and Pb = ^ and get 
the Hamiltonian H — apa + hpi, — L. One of its canonical equations is a = 6 
and another one is equivalent to the fourth-order Euler-Lagrange equation 
following from L. By these properties, L^^y is uniquely determined up to the 
integration constant. Contrary to other approaches, no constraint is needed. 



3 Fourth— order gravity 

In Rainich (1925, ref. [17]) the electromagnetic field was calculated from the 
curvature tensor. This was cited in Kuchaf (1963, ref. [18]) as example for 
the geometrization programme; in [18] on meson fields (now called scalar 
fields), Kuchaf gives a kind of geometrization by using a relation between '0 
and i?, then he gets the equation 

OR R=0 

2 

which is of fourth order in the metric. It looks like fourth-order gravity as 
we are dealt with, but he does not deduce it from a curvature squared action. 
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3.1 Duality theorems 

In Bekenstein (1974, ref. [19]) the conformal transformation from Einstein's 
theory with a minimally coupled (0) to a conformally coupled {ip) scalar field 
is proven where additional conformally invariant 

matter (radiation) is allowed. For SnG — 1 one has 

V' = v^tanh(0/V6) 

If radiation is absent then it works also with "coth" instead of "tanh". This 
is reformulated in his theorem 2: If Qij and ip form an Einstein-conformal 
scalar solution, then Qij = ^ip'^Qij and -ijj = form a second one. One can 
see that this is a dual map because by applying the operator " twice, the 
original solution is re-obtained. 

Let us comment this theorem 2: For the conformal scalar field one has 
the effective gravitational constant Ge// defined by 

1 = I _t. 
SnGeff 6 

A positive value G^ff implies a negative value Geff- By changing the overall 
sign of the Lagrangian one can achieve a positive effective gravitational con- 
stant at the price of the scalar field becoming a ghost (wrong sign in front 
of the kinetic term). So, Bekenstein has given a conformal transformation 
from Einstein's theory with a conformally coupled ordinary scalar field to 
Einstein's theory 

with a conformally coupled ghost. From this property one can see that 
this duality relation is different from the duality theorem to be deduced at 
the end of this subsection, because there one has effectively ordinary scalar 
fields at both sides. 

Later but independently of [19] the conformal equivalence between min- 
imally and conformahy coupled scalar fields with Ge// > was generalized 
in [20] by the inclusion of several self-interaction terms. 

The conformal transformation from fourth-order gravity to Einstein's 
theory with a minimally coupled scalar field was deduced in several steps: 
Bicknell (1974, ref.[21]) found the transformation for L = E?; the conformal 
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factor is R, and after the transformation one gets Einstein's theory with non- 
vanishing A-term and a massless minimally coupled scalar field as source. 
Next steps sec e.g. [22]. In [23] besides the conformal transformation it is 
shown that the trace ~ leads to a term like R^ In R in L. Ref. [24] 
generalizes it 

by the inclusion of non-minimally coupled scalar fields. Jakubicc and 
Kijowski (1988/89, ref. [25]) generalize the conformal transformations to 
more general transformations of the metric. 

Buchdahl (1978, ref. [26]) showed: For L = R^ the conformal factor R^ 
(if it is 7^ 0) transforms solutions to solutions and represents a dual map in 
the set of solutions. (Another conformal factor than in [21]!) [27] generalizes 
this dual map to other non-linear Lagrangians L{R), the conformal factor 
being (^)^- (Again, this conformal factor is the square of that conformal 
factor which is necessary to transform to Einstein's theory with a minimally 
coupled scalar field.) 

A further type of transformations was presented by Buchdahl in 1959, 
ref. [28]. For a space-time T4 possessing a non-null hyp ersurf ace-orthogonal 
Killing vector one can produce solutions with a scalar field as follows. Let 
the Killing vector be the metric is g^i with g^i^i = 0, gu ^ and g\oi — 
where greek indices take all values except 1. Let R^i — and the dimension 
n > 2 but n 7^ 3. We fix two reals A and B and define a new metric g^ by 
~9\o. = 0, 

(We consider the case gu > 0, the other sign is treated analogously.) This 
is a conformal transformation for A = B + 1. Defining gu = e^^ one has 
= iff (= if and only if) 0-0 = 0. So, □'0 = is presumed from the 
beginning. 

Further, it holds: ~ □■0 = is then identically fulfilled iff 

A = 1 - S(n - 3) 

(For n = 4, this means A = 1 — B.) 

B = 0, A = 1 represents the identical transformation. For non-identical 
transformations the equation is compatible with a conformal transformation 
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for n = 2 only. This is in agreement with the fact that the D'Alembert 
operator is conformally invariant for n = 2 only. 

Buchdahl's results are: For n — 4 and A = —1, i.e., B = 2, one gets 
a vacuum solution which is in general different from the initial one. For 
n — 4 and 1^41 ^ 1, one gets a solution of Einstein's theory with a massless 
minimally coupled scalar field (which is proportional to ip) as source. This 
supplements [21], because here no A-term is needed. 

In [29, 30] the problem is discussed which of the conformally equivalent 
metrics in these theorems is the physical metric. Magnano and Sokolowski 
(1994, ref. [30]) represents a good review to this theme. 

Let us now deduce the duality theorem announced in the introduction 
which shall close a gap in the set of the 
aforementioned results. 

Let 

9ij = Qij (3.1) 

and Dc = n - |. Then 

= (3.2) 

reflects the conformal invariance of the operator O^. if applied to a scalar. 
The following consideration is restricted to the case that both R and R are 
positive; the other sign can be dealt analogously. 

We put U = InR into eq. (3.1) and apply the identity (3.2) to the 
constant scalar = 1. We multiply by (— 6i?) and get the identity 

RR ^ 1 - —DR (3.3) 
R"^ 

Further, we define the operator " given by 

9ij = 9ij (3.4) 

to be a dual one if it coincides with its inverse operator. Then the following 
conditions are equivalent: 
1. " is a dual operator. 
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2. RR ^ 1 

3. OR^O 

4. ^0^ = 

Proof. 1. -^2. is simply the explicit form of the definition of duahty. 
2. 44> 3. follows from eq. (3.3). 3. 44> 4. is a consequence of the duality 
property. 

Remcirks. 1. The conformal factor R^ in eq. (3.4) is crucial for the 
validity of the duality. Already for spaces with constant curvature scalar one 
can see: R—^ requires the conformal factor to be R^. 2. This consideration 
was inspired by Buchdahl's paper [26] from 1978. 

Now we are prepared to show a duality relation between pairs of fourth- 
order gravity theories. (We formulate it only for dimension 4, other dimen- 
sions > 2 give similar resuhs.) Let L = F{R) with G = '^y^Q, H=^=^0 
and gij — gij. If gij is a solution of the fourth-order equation following 
from L then it holds 

R = 3R/G^ - AF/G^ (3.5) 
One has ^ 7^ iff (= if and only if) 

^ 6H{2F-GR) (3.6) 

Remark. If i? is a constant then 2F = GR follows from the field equation, 
(3.6) is automatically fulfilled, and (3.5) reduces to the known relation R = 
R/G\ 

Let (3.6) be fulfilled in the following. Then eq. (3.5) can be inverted 
locally as i? = R{R). We define 

F = 2R/G^ - 3F/G^ (3.7) 

where R = R{R) has to be inserted into the r.h.s. 

The following theorem holds under the presumptions formulated above. 

Theorem. L — F{R) defines a fourth-order theory dual to L = 
F{R) and Qij is a solution of its field equation. 

Duality means that by applying this procedure twice the original theory 
with original solution gij is obtained. 
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Remcirks. 1. The most problematic step in formulating this theorem 
was to find the formulas (3.5) and (3.7); the existence of such a theorem was 
announced in [27], but these two crucial formulas are presented here for the 
first time. 

2. We call two theories F, F to be similar (i.e., they go into each other 
by a change of the length unit) if there exist non-vanishing reals a and fi 
such that F{R) — aF{j3R). It holds: Similar theories have similar duals. 

Proof of the Theorem. The trace of the field equation following from 
L reads 

GR +3DG = 2F 

It must be supplemented by the condition that the trace-free part of the 
tensor GRij — G-ij vanishes. We use identities analogous to eq. (3.3) with G 
instead of R and the relation holding for the r.h. sides of eqs. (3.5) and (3.7): 

This ensures the validity of 

G 

The rest of the proof is lengthy but straightforward. 

Examples. 1. Let L — ^ with m 7^ 0, 1. Unequality (3.6) excludes 
m = I and m = |. For all remaining reals m one gets the dual L — cR^ 
with 771 = and a constant c = c(m), cf. [27]. 

2. Let L — R — aR^. The inversion of eq. (3.5) is not possible in closed 
form, but in the vicinity of fiat space we may expand into powers of R and 
get up to third order 

L^R- ak^- Aa^R^ 

Result: At least a special type of cubic terms in the Lagrangian can be 
absorbed by a suitable conformal tranformation. 



dR _ dF 
dR ~ dR 



dR G 
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3.2 Instability of i?^- theories 

This subsection deals with the classical instabihty of fourth-order theories 
following from a non-linear Lagrangian L[R). 

(Quantum instabilities will be commented in set. 5.) 

Teyssandier and Tourrenc (1983, ref. [31]) solved the Cauchy-problem 
for this theory, let us shortly repeat the main ingredients. 

The Cauchy problem is well-posed (a property which is usually required 
to take place for a physically sensible theory) in each interval of i?-values 
where both dL/dR and d'^L/dF^ are different from zero. The constraint 
equations are similar as in General Relativity: the four Oi-component equa- 
tions. What is different are the necessary initial data to make the dynamics 
unique. More exactly: Besides the data of General Relativity one has to 
prescribe the values of R and ^ at the initial hypersurface. This coin- 
cides whith the general experience: Initial data have to be prescribed till the 
highest-but-one temporal derivative appearing in the field equation (here: 
fourth-order field equation, ^ contains third-order temporal derivatives of 
the metric). Under this point of view, classical stability of the field equation 
means that a small change of the Cauchy data implies also a small change of 
the solution. 

Now we are prepared to classify the stability claims found in refs. [32 
- 36]. To simplify we specialize to the Lagrangian L = R - eR^ with the 
non-tachyonic sign e > and restrict to the range ^ > 0, i.e. R < ^ 

On the one hand, refs. [32, 33] find a classical instability of the Minkowski 
space-time for this case. (Mazzitelli and Rodrigucs [33] cite Gross, Perry and 
Yaffe (1982, ref. [34]) with the sentence "The Minkowski solution in general 
relativity has been proven to be stable." which refers to the positive energy 
theorem of general relativity.) 

On the other hand, refs. [35, 36] find out that the Minkowski space-time 
is here not more unstable than in General Relativity itself. What looks like 
a contradiction from the first glance is only a notational ambivalence as can 
be seen now: The main argument in refs. [32, 33] is that an arbitrarily 
large value ^ is compatible with small values of if^ and R^. In refs. [35, 
36] however, following the Gauchy-data argument [31], (^ being part of 
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the Cauchy data which are presumed to be small) stabihty of the Minkowski 
space-time is obtained in the version: If the Cauchy data are small (meaning: 
close to the Cauchy data of the Minkowski space-time) then the fourth-order 
field equation bounds the solution to remain close to the Minkowski space- 
time. 

The argument of ref. [35] is a little bit different: There the conformal 
transformation to Einstein's theory with a scalar field $ [22] is applied; it 
is observed that in the F(i?) -theory there are never ghosts which implies 
stability. Now, $ and ^ belong to the Cauchy data which is equivalent to 
the data R, ^ in the conformal picture thus supporting the Cauchy data 
argument given at the beginning of this subsection. 

4 Cosmology 

Several papers [37 - 43] apply the conformal transformation theorem [21, 22] 
to cosmology; so for interpreting the cosmological singularity [37], for dealing 
with anisotropic models [38], with transformation to Brans-Dicke extended 
inflation [40]. Ref. [42] mentions that < 1 is possible even if A; = 1. The 
other ones apply the theorem mainly as a mathematical device to transform 
solutions to solutions of the other theory. 

In the following two subsections we discuss exact solutions directly found 
for fourth-order gravity (subsection 4.1: a spatially flat Priedmann model, 
subsection 4.2: a Kantowski-Sachs model). 

4.1 Friedmann models 

In 1988, Chimento [44] found an exact solution for fourth-order gravity in 
a spatially flat Priedmann model. He also found out that in the tachyonic- 
free case the asymptotic matter-dominated Priedmann solution is stable, and 
no fine-tuning of initial conditions is necessary to get the final (oscillating) 
Priedmann stage; particle production of non-conformal fields may backreact 
to damp the oscillations. 
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[45] generalizes [44]: here the Dirac equation is considered, the result 
is that there appear spinor field oscillations and the qualitative behaviour 
remains essentially the same. 

Let us present the exact solution of [44] . For the spatially flat Friedmann 
model with Hubble parameter H = a/ a he solves the fourth-order field 
equation with vacuum polarization term. The zero-zero component equation 
reads 

2HH - + m^H + -H^ + H'^ ^0 (4.1) 

4 

The if^-term stems from the vacuum polarization and the H^-term from the 
Einstein tensor. The remaining ingredients of eq. (4.1) come from the term 

in the Lagrangian. (Here we only present the tachyonic-free case with 
A = and | in front of if^.) The factor in front of if^ should not infiuence 
the weak-field behaviour because for if this factor only changes the 
effective gravitational constant. 

From eq. (4.1) the discussion of subsection 3.2 becomes obvious: (4.1) 
represents a third-order equation for the cosmic scale factor a; it is a con- 
straint and not a dynamical equation. (It is only due to the high symmetry, 
that accidentally the validity of the constraint implies the validity of the dy- 
namical equation.) Supposed, eq. (4.1) would be the true dynamical equation 
for a theory, then the instability argument of [32] could apply. 

The ansatz for solving eq. (4.1) 

leads to a non-linear third-order equation for s 

2ss^^') -s^ + s"^ = (4.2) 

Derivative with respect to t yields the equation s^^^ -|- s = being linear in s 
and having the solution 

s = ci + C2t + C3 sin(t + C4) 

Inserting this solution into the original equation gives the restriction |c2| = 
|c3|. Let us discuss this solution: C2 — leads to the uninteresting fiat space- 
time. So, now let C2 7^ 0. Adding tt to C4 can be absorbed by a change of 
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the sign of C3. Therefore, C2 = C3 without loss of generahty. Multiphcation 
of s by a constant factor does not change the geometry, so let C2 = 1. A 
suitable time-translation leads to Ci = 0. Finally, the cosmic scale factor is 
calculated as a = s^l^ leading to 

2 

a = [t + sin(t + 04)]^/^ - t^/^ [1 + — sin(t + C4)] (4.3) 

The r.h.s. of eq. (4.3) gives in an elegant way the late-time behaviour 
already deduced in [46]. The factor 1/t in front of the "sin" -term shows that 
the oscillations due to the higher-order terms are damped. The total energy 
"sitting" in these oscillations, however, remains constant in time (because of 
the volume-expansion) cf. Suen (1994, ref. [32]) and can be converted into 
classical matter by particle creation. 

Let us mention some further cosmological solutions with higher-order 
gravity: [47] discusses the L(i?)-stabihty with a conformally coupled scalar 
field. Ref. [48] (partial results of it can be found in [49]) deals with fourth- 
order cosmological models of Bianchi-type I and power-law metrics, i.e. 

ds^ ^dt^-j^ {dx'f 

i=l 

with real parameters pj. The suitable notation 

3 

i=l 

gives the following: ai = a2 = 1 is the usual Kasner solution for Einstein's 
theory, + 02 = 2ai is the condition to be fulfilled for a solution in L — R'^. 
Refs. [50, 51] also 

discuss i?^-models. Briining, Coule and Xu (1994, ref. [36]) consider 
inflationary cosmology with a Lagrangian 

L ^ R+ XRi,^Rf"'/R 

and mention that it is unclear under which circumstances the existence of 
the Weyl term in anisotropic models allows the de Sitter space-time to be an 
attractor solution. Ref. [52] deals with anisotropic Bianchi-type IX solutions 
for L — R^. They look for chaotic behaviour analogous to the mixmaster 
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model in Einstein's theory. Ref. [53] gives exact solutions for L — B? and a 
closed Priedmann model, ref. [54] discusses the bounce in closed Priedmann 
models for L = R — eB?. Supplementing the discussion of [54, eq.(l)] let 
us mention: In the non-tachyonic case, there exist periodically oscillating 
models with an always positive scale factor a. Ref. [55] looks for chaos in 
isotropic models, e.g. by conformally coupled massive scalar fields in the 
closed universe. The papers [56, 57] consider the stability of power-law 
inflation for L = R™- within the set of spatially fiat Friedmann models. Refs. 
[58] give overviews on higher-order cosmology, especially chaotic inflation 
as an attractor solution in initial-condition space. [59] deals with quantum 
gravitational effects in the de Sitter space-time, and [60] gives a classification 
of inflationary Einstein-scalar-field-models via catastrophe theory. Ref. [61] 
considers Chern-Simon terms in Bianchi cosmologies and the cosmic no-hair 
conjecture. The axion with field strength Hijk puts an extra hair on black 
holes. Its energy-momentum tensor does not fulfil the energy conditions, and 
so one gets both recoUapsing solutions and ever-expanding solutions which 
are essentially anisotropic also for late times. 

4.2 Kantowski— Sachs models 

Before we come to the fourth-order solution by Accioly let us mention some 
results on Kantowski-Sachs models in general. 

The solution found in 1950, cf. ref. [62] - now it is called Nariai solution 
- is the only static spherically symmetric solution of the Einstein equation 
with positive A-term which cannot be written in Schwarzschild coordinates. 
It has a six-dimensional isometry group and is of Kantowski-Sachs type; it 
represents the direct product of two two-dimensional spaces of equal and 
non-vanishing constant curvature (in short: S"^ x S'^). One of the many 
possibilities to present it is 

ds^ - (1 - Kr^)dt' - YZrj^ ~ ('^^^ + (^-^^ 

Ref. [63] discusses the stability of the Bertotti-Robinson (also the direct 
product of two two-dimensional spaces of non-vanishing constant curvature. 
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but with vanishing 4-curvature scalar, in short S"^ x S'^) and Nariai solu- 
tions. Ref. [64] is an obituary to H. Nariai, it especially mentions the Nariai 
solution. Kofman, Sahni and Starobinsky (1983, ref. [65]) get the result 
that there is no particle production in the Nariai solution. [66] discusses 
the analogous question for the Bertotti-Robinson metric, the authors of ref. 
[67] consider the Nariai metric and its decay into de Sitter and Kasner-like 
space-times. They consider essentially Einstein's vacuum theory with posi- 
tive A-term and Kantowski-Sachs metric. 

Torrcnce and Couch (1988, ref. [68]) showed how the de Sitter space-time 
can be presented locally as a Kantowski-Sachs cosmological model. More- 
over, no other Robertson- Walker space-time can be presented as Kantowski- 
Sachs model, cf. also [69] for this question. 

Moniz (1993, ref. [70]) discussed the Kantowski-Sachs models in Ein- 
stein's theory with positive A-term in relation to the no-hair conjecture and 
got the following result. The majority of solutions is asymptotically de Sit- 
ter, a small number recoUapses: infinite to finite measure if dB dB is taken 
as measure in the initial condition space, and B is the radius of the S'^ of the 
model: 

ds^ = dt" - A\t)dr' - B\t){de^ + sin' ^#') 

Accioly (1987/88, rcfs. [71, 72]) found an exact cosmological solutions of the 
Godel type for higher-order gravity, the new solution he found is the direct 
product of a 3-space of non-vanishing constant curvature with the real line in 
space-like direction (see [72, eq. (18)]). This gives a 7-dimensional isometry 
group. 

Comment: Because of the high symmetry, this exact solution belongs 
both to the Godel and to the Kantowski-Sachs classes of metrics (with A — 
const, and B — coshi). The latter class is more popular, so we classify it 
here and not primarily as Godel model. 

The direct product of a 3-space of positive constant curvature with the 
real line in time-like direction is known as Einstein's static universe. So, 
by imaginary coordinate transformations, both metrics can be mapped onto 
each other locally. This imphes that Accioly's solution is conformally flat 
(because Einstein's universe carries this property). In the set of conformally 
flat space-times, the term RijE}^ gives the same variational derivative as |i?'. 
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so Accioly's paper has to be classified in the set of non-hnear Lagrangians 
L — F{R); he assumes F to be a quadratic function of i?. As one knows 
there is a critical value of the curvature scalar in such theories, it is defined 
by ^ = 0. (Cf. subsection 3.2 above: there the Cauchy problem is not 
well-posed.) 

At these values, the fourth-order difi^erential equation reduces to the sec- 
ond order one F = 0; it turns out that Accioly's solution obeys this critical 
value so it is unstable similar to Einstein's static universe. More directly: Let 
L = [R — RoY with a constant Ro, then each solution of the second-order 
equation R = Rq solves also the fourth-order field equation following from 
L. So a lot of solutions (including Accioly's one) can be found, but they all 
live in the region where the Cauchy problem is ill-posed. 

5 Summary 

The scope of this paper was to present the foundations necessary to deduce 
the Wheeler-de Witt equation for a cosmological minisuperspace model in 

fourth-order gravity. 

The method (sporadically developed in [51] for L = R? and a spatially 
flat Friedmann model) to handle with eqs. (1.3 - 1.6) was systematically 
generalized in set. 2 to give a Hamiltonian formulation of a general fourth- 
order theory. The possibility of deducing this method makes it clear that 
the method of ref. [51] is not restricted to highly symmetric models. The 
alternate Hamiltonian formulation has some advantages in comparison with 
Ostrogradski's one: No constraint is needed, the Hamiltonian is typically a 
quadratic function in the momenta. (Ostrogradski's approach leads always 
to a Hamiltonian linear in the momenta which gives artificial factors i in the 
Schrodinger equation.) The calculation of the momenta from the Lagrangian 
follows the usual equations (2.8) whereas the Ostrogradski approach needs 
some additional terms. Our approach is less ambiguous, cf. eq. (2.7). 

One could pose the question whether both approaches are equivalent on 
another level, this is not fully excluded, but even if it is the case, the approach 
deduced here is more directly applicable to fourth-order quantum cosmology. 
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The fact that the Jacobian eq. (2.28) is always negative excludes the 
possibility to get a positive definite Jacobi metric in eq. (1.2). This is 
one of the many possibilities to say what is meant by the phrase "fourth- 
order theories are always unstable". The Jacobi metric plays the role of 
the conformally transformed superspace-metric used in quantum cosmology. 
And here the circle can be closed: In Einstein's theory (both for Lorentzian 
and Euclidean signature of the underlying manifold) the superspace-metric 
has Lorentzian signature and cannot be positive definite. So we get once more 
the result of subsection 3.2: Fourth-order gravity contains some instabilities, 
but only those which it has in common with General Relativity. 

To decide the quantum instability of the Minkowski or de Sitter space- 
times in fourth-order gravity one must solve the corresponding Wheeler-de 
Witt equations (Mazzitelli and Rodrigues [33] deduced them for the spatially 
fiat Friedmann model and the Lagrangian L = R — eR?) and has to interpret 
them carefully. This has to be done yet. 

In [50] it is mentioned that a classical theory with higher derivatives has 
instabilities: "At the quantum level, the difference is even more dramatic. 
Noncommuting variables in the lower-derivative theory, such as position and 
velocities, become commuting in the higher-derivative theory." Remark of 
U. Kasper to this sentence: "The uncertainty relation is primarily between 
positions and momenta. If the momentum is independent of the velocity 
then commuting position and velocity need not bother." 

The duality theorem deduced in subsection 3.1 is a method to construct 
new solutions of fourth-order gravity from known solutions of a (possibly 
other) fourth-order theory. It gives non-trivial results only for solutions with 
a non-constant curvature scalar, e.g. the Chimento solution [44] which has 
been rewritten in subsection 4.1. Kantowski-Sachs models (whose quantum 
cosmology is considered in [73]) have begun to discuss in subsection 4.2; this 
shall be completed elsewhere. 
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